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The STAR Collaboration at RHIC has measured two-pion correlation functions from p + p collisions at 
= 200 GeV . Spatial scales are extracted via a femtoscopic analysis of the correlations, though this analysis 
is complicated by the presence of strong non-femtoscopic effects. Our results are put into the context of the 
world dataset of femtoscopy in hadron-hadron collisions. We present the first direct comparison of femtoscopy 
in p + p and heavy ion collisions, under identical analysis and detector conditions. 



I. INTRODUCTION AND MOTIVATION Quantum Chromodynamics (QCD) from numerous direc- 

tions. The extraordinary flexibility of the machine permits 
collisions between heavy and light ions at record energies (up 
The experimental program of the Relativistic Heavy Ion to = 200 GeV ), polarized and unpolarized protons, and 
ColUder (RHIC) at Brookhaven National Laboratory probes strongly asymmetric systems such as c/ + Au . The proton col- 
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lisions are the focus of an intense program exploring the spin 
structure of the nucleon. However, these collisions also serve 
as a critical "baseline" measurement for the heavy ion physics 
program that drove the construction of RHIC. 

Studies of ultrarelativistic heavy ion collisions aim to ex- 
plore the equation of state of strongly interacting matter. The 
highly dynamic nature of the collisions, however, does not 
allow a purely statistical study of static matter as one might 
perform in condensed matter physics, but rather requires a de- 
tailed understanding of the dynamics itself. If a bulk, self- 
interacting system is formed (something that should not be as- 
sumed a priori), the equation of state then plays the dynamic 
role of generating pressure gradients that drive the collective 
expansion of the system. Copious evidence d-Q] indicates 
that a self-interacting system is, in fact, generated in these col- 
lisions. The dynamics of the bulk medium is reflected in the 
transverse momentum (pr) distribution 15ijfi| and momentum- 
space anisotropy (e.g. "elliptic flow") Qlsll of identified par- 
ticles at low pT- These observables are well-described in a 
hydrodynamic scenario, in which a nearly perfect (i.e. very 
low viscocity) fluid expands explosively under the action of 
pressure gradients induced by the collision i^. 

Two-particle femtoscopy ifToll (often called "HBT" anal- 
ysis) measures the space-time substructure of the emitting 
source at "freeze-out," the point at which particles decouple 
from the system [e.g. [HI]. Femtoscopic measurements play 
a special role in understanding bulk dynamics in heavy ion 
collisions, for several reasons. Firstly, collective flow gener- 
ates characteristic space-momentum patterns at freezeout that 
are revealed ifTTIl in the momentum-dependence of pion "HBT 
radii" (discussed below), the transverse mass dependence of 
homo gen eity lengths III2I1 . and non-identical particle correla- 
tions lUOl VIM . Secondly, while a simultaneous description 
of particle-identified pr distributions, elliptic flow and femto- 
scopic measurements is easily achieved in flow-dominated toy 
models [e.g.i^, achieving the same level of agreement in a re- 
alistic transport calculation is considerably more challenging. 
In particular, addressing this "HBT puzzle" fl4}l has led to a 
deeper understanding of the freezeout hypersurface, collectiv- 
ity in the initial stage, and the equation of state. Femtoscopic 
signals of long dynamical timescales expected for a system 
undergoing a first-order phase transition ifisi IT6I1 . have not 
been observed ifTlll . providing early evidence that the system 
at RHIC evolves from QGP to hadron gas via a crossover IitIi . 
This sensitive and unique connection to important underlying 
physics has motivated a huge systematic study of femtoscopic 
measurements in heavy ion collisions over the past quarter 
century ifTHl . 

HBT correlations from hadron (e.g. p + p) and lepton (e.g. 
e+ + ) collisions have been extensively studied in the high 
energy physics community, as well |[l8l - t2C)ll . although the the- 
oretical interpretation of the results is less clear and not well 
developed. Until now, it has been impossible to quantitatively 
compare femtoscopic results from hadron-hadron collisions 
to those from heavy ion collisions, due to divergent and often 
undocumented analysis techniques, detector acceptances and 
fitting functions historically used in the high energy commu- 
nity 1I20I1 . 



In this paper, we exploit the unique opportunity offered by 
the STAR/RHIC experiment, to make the first direct com- 
parison and quantitative connection between femtoscopy in 
proton-proton and heavy ion collisions. Systematic compli- 
cations in comparing these collisions are greatly reduced by 
using an identical detector and reconstruction software, colh- 
sion energies, and analysis techniques (e.g. event mixing ll2lll . 
see below). We observe and discuss the importance of non- 
femtoscopic correlations in the analysis of small systems, and 
put our femtoscopic results for p + p collisions into the con- 
text both of heavy ion collisions and (as much as possible) 
of previous high-energy measurements on hadron-hadron and 
e++e^ collisions. These results may play a role in un- 
derstanding the physics behind the space-momentum corre- 
lations in these collisions, in the same way that comparison 
of p + p and heavy ion collision results in the high-pj- sec- 
tor is crucial for understanding the physics of partonic energy 
loss lfll-l4l l22ll . Our direct comparison also serves as a model 
and baseline for similar comparisons soon to be possible at 
higher energies at the Large Hadron Collider. 

The paper is organized as follows. In Sectionllll we discuss 
the construction of the correlation function and the forms used 
to parameterize it. Section|III]discusses details of the analysis, 
and the results are presented in Section HV] In Section |V] we 
put these results in the context of previous measurements in 
Au + Au and p + p{p) collisions. We discuss the similarity 
between the sy stematics of HBT radii in heavy ion and particle 
collisions in Section[yI]and summarize in Section lVTll 

II. TWO-PARTICLE CORRELATION FUNCTION 

The two-particle correlation function is generally defined as 
the ratio of the probability of the simultaneous measurement 
of two particles with momenta pi and p2, to the product of 
single-particle probabilities. 



Cipi^pi) ■ 



PjPUPl) 
P{PX)P{P2)' 



In practice, one usually studies the quantity 



Cp {q) 



Ap jq) 
BpiqY 



(1) 



(2) 



where q = pi — p2 is, the relative momentum. A{q) is the 
distribution of the pairs from the same event, and B{q) is 
the reference (or "background") distribution. B contains all 
single-particle effects, including detector acceptance and ef- 
ficiency, and is usually calculated with an event-mixing tech- 
nique llllll2lll . The explicit label P (= {pi + p2) /2) empha- 
sizes that separate correlation functions are constructed and 
fitted (see below) as a function of q, for different selections 
of the total momentum P; following convention, we drop the 
explicit subscript below. Sometimes the measured ratio is nor- 
malized to unity at large values of \q\; we include the normal- 
ization in the fit. 

In older or statistics-challenged experiments, the cor- 
relation function is sometimes constructed in the o ne- 

dimensional quantity Qi^v = \J {pi — P2)^ — {E\ — El)" or 
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two-dimensional variants (see below). More commonly in re- 
cent experiments, it is constructed in three dimensions in the 
so-called the "out-side-long" coordinate system lEsl - Esll . In 
this system, the "out" direction is that of the pair transverse 
momentum, the "long" direction is parallel to the beam, and 
the "side" direction is orthogonal to these two. We will use 
the subscripts "o," "/" and "i" to indicate quantities in these 
directions. 

It has been suggested ll26l - l28ll to construct the three- 
dimensional correlation function using spherical coordinates 

^„ = 1^1 sin0cos(|), = 1^1 sin 9 sin (|), ^/ = |^|cos9. 

(3) 

This aids in making a direct comparison to the spatial sepa- 
ration distribution through imaging techniques and provides 
an efficient way to visuaUze the full three-dimensional struc- 
ture of C{q). The more traditional "Cartesian projections" 
in the "o," "s" and "/" directions integrate over most of the 
three-dimensional structure, especially at large relative mo- 
mentum ifTilEsli. 

Below, we will present data in the form of the spherical 
harmonic decomposition coefficients, which depend explicitly 

on 1^1 as 



Ai,A\q\) 



1 



d^d{co& Q)C{\q\,Q,^) Yijn (9, (^) . (4) 



The coefficient Aoo(|g'|) represents the overall angle- 
integrated strength of the correlation. A20 (|^|) and A22 {\q\) 
are the quadrupole moments of C at a particular value of |^| . In 
particular, A22 quantifies the second-order oscillation around 
the "long" direction; in the simplest HBT analysis, this term 
reflects non-identical values of the /?„ and Rs HBT radii (c.f. 
below). Coefficients with odd I represent a dipole moment 
of the correlation function and correspond to a "shift" in the 
average position of the first particle in a pair, relative to the 
second 126142811 . In the present case of identical particles, the 
labels "first" and "second" become meaningless, and odd-/ 
terms vanish by symmetry. Likewise, for the present case, 
odd-m terms, and all imaginary components vanish as well. 
See Appendix B of [|28il for a full discussion of symmetries. 

In heavy ion collisions, it is usually assumed that all of the 
correlations between identical pions at low relative momen- 
tum are due to femtoscopic effects, i.e. quantum statistics and 
final-state interactions llllll . At large |^|, femtoscopic effects 
vanish [e.g. |Tl]]. Thus, in the absence of other correlations, 
C (q) must approach a constant value independent of the mag- 
nitude and direction of q; equivalently. A/ ,„ (|^|) must vanish 
at large \q\ for I 7^ 0. 

However, in elementary particle collisions additional struc- 
ture at large relative momentum (|^| > 400 MeV/c) has been 
observed [e.g.i20, 29- 33.1 . Usually this structure is parameter- 
ized in terms of a function Q. {q) that contributes in addition to 
the femtoscopic component Cp (q). Explicitly including the 
normalization parameter , then, we will fit our measured 
correlation functions with the form 



C{q)^:^^■CF{q)■^iq). 



(5) 



Below, we discuss separately various parameterizations of the 
femtoscopic and non-femtoscopic components, which we use 



in order to connect with previous measurements. A historical 
discussion of these forms may be found in ll20ll . 

We use a maximum-likelihood fit to the correlation func- 
tions, though chi-square minimization yields almost identical 
results, and we give the values for all fits below. As we 
shall see, none of the functional forms perfectly fits the data. 
However, the characteristic scales of the source can be ex- 
tracted and compared with identical fits to previous data. 

A. Femtoscopic correlations 

Femtoscopic correlations between identical pions are dom- 
inated by Bose-Einstein symmetrization and Coulomb final 
state effects in the two-pion wavefunction ifTHl . 

In all parameterizations, the overall strength of the fem- 
toscopic correlation is characterized by a parameter A, ifTTIl . 
Historically called the "chaoticity" parameter, it generally ac- 
counts for particle identification efficiency, long-Uved decays, 
and long-range tails in the separation distribution ll34ll . 

In the simplest case, the Bose-Einstein correlations are of- 
ten parameterized by a Gaussian, 



(6) 



Kopylov and Podgoretskii 113511 introduced an alternative, 
two-dimensional parameterization 



where R\nv is a one dimensional "HBT radius." 



CpiqT.qo) = 1 +A. 



2/1 (qiRB) 



qjRB 



(7) 



where qj is the component of q orthogonal to P, qo = Ei — 
E2, Rb and 1 are the size and decay constants of a spherical 
emitting source, and J\ is the first order Bessel function. This 
is similar to another common historical parameterization [e.g. 
[36I1 characterizing the source with a spatial and temporal scale 



Cpiq^qo) 



1 



(8) 



Simple numerical studies show that Rq from Eq. [Sjis ap- 
proximately half as large as, Rb obtained from Eq. 171 1I20I l3a. 

With sufficient statistics, a three-dimensional correlation 
function may be measured. We calculate the relative mo- 
mentum in the longitudinally co-moving system (LCMS), in 
which the total longitudinal momentum of the pair, pi i + 
p/,2, vanishes Issll . For heavy ion and hadron-hadron colli- 
sions, this "longitudinal" direction I is taken to be the beam 
axis ifTlll : for e+ + collisions, the thrust axis is used. 

For a Gaussian emission source, femtoscopic correlations 
due only to Bose-Einstein symmetrization are given by [e.g. 



CAqo,qs,qi) = l+y^e-^"^'-^'''"-^>'''l 



(9) 



where Rg, Rs and Ri are the spatial scales of the source. 

While older papers sometimes ignored the Coulomb final- 
state interaction between the charged pions ll20ll . it is usually 
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included by using the Bowler- Sinyukov 11391 14011 functional 
form 

Cf (e,nv) = (1 - ^) + X/:coul (Ginv) (l +e-eL«L) , (10) 

and in 3D, 

CF{qo,qs,qi) = (1 - A.) + A-^Tcoul (Ginv) 

X (^1 +e-'?"*»-'i'?^?-''r^^^ . (11) 

Here, K^oui is the squared Coulomb wavefunction integrated 
over the source emission points and over the angles of the 
relative momentum vector in the pair rest frame. 

B. Non-femtoscopic correlations 

In the absence of non-femtoscopic effects, one of the forms 
for Cf (q) from Section UTAI is fitted to the measured correla- 
tion function; i.e. H = 1 in Equation|5] Such a "standard fit" 
works well in the high-multiplicity environment of heavy ion 
collisions fll]]. In hadron-hadron or e+ + collisions, how- 
ever, it does not describe the measured correlation function 
well, especially as |<7| increases. Most authors attribute the 
non-femtoscopic structure to momentum conservation effects 
in these small systems. While this large-|^| behavior is some- 
times simply ignored, it is usually included in the fit either 
through ad-hoc [30] or physically-motivated |f28l l terms. 

In this paper, we will use three selected parameterizations 
of the non-femtoscopic correlations and study their effects on 
the femtoscopic parameters obtained from the fit to experi- 
mental correlation functions. The first formula assumes that 
the non-femtoscopic contribution can be parameterized by a 
first-order polynomial in ^-components (used e.g. in ll4ll - l45ll ). 
Respectively, the one- and three-dimensional forms used in 
the literature are 



£2(^) = l+5^ 



(12) 



and 



Q.{q) = ^{qc,q.„qi) = I + 5oqo + 5,qs + 8iqi. (13) 



For simplicity, we will use the name "8 — q fit" when we fit 
Eq. [12] or [13] to one- or three-dimensional correlation func- 
tions. 

Another form ll46ll assumes that non-femtoscopic correla- 
tions contribute |^| -independent values to the / = 2 moments 
in EquationlH In terms of the fitting parameters ^ and p, 

n ( 1^1, cos e, (j)) = i2 (cose, (j)) = 

l+2\^(pF2,o (cose,(j))+2CRe[y2,2 (cose,(j))]) = 




(3cos2e-l) + ^ 



^^sin2ecos2(^. (14) 



For simplicity, fits using this form for the non-femtoscopic 
effects will be referred to as — p fits." 

These two forms (as well as others that can be found in 
literature ll20ll ) are purely empirical, motivated essentially by 



the shape of the observed correlation function itself. While 
most authors attribute these effects primarily to momentum 
conservation in these low-multiplicity systems, the parame- 
ters and functional forms themselves cannot be directly con- 
nected to this or any physical mechanism. One may iden- 
tify two dangers of using an ad-hoc form to quantify non- 
femtoscopic contributions to C{q). Firstly, while they de- 
scribe (by construction) the correlation function well at large 
\q\, for which femtoscopic contributions vanish, there is no 
way to constrain their behaviour at low |^| where both femto- 
scopic and (presumably) non-femtoscopic correlations exist. 
Even simple effects like momentum conservation give rise to 
non-femtoscopic correlations that vary non-trivially even at 
low 1^1 . Misrepresenting the non-femtoscopic contribution 
in D. [q) can therefore distort the femtoscopic radius param- 
eters in Cf (q), especially considering the small radius val- 
ues in p + p collisions. Secondly, there is no way to estimate 
whether the best-fit parameter values in an ad-hoc functional 
form are physically "reasonable." 

If the non-femtoscopic correlations are in fact dominated by 
energy and momentum conservation, as is usually supposed, 
one may derive an analytic functional form for D.. In partic- 
ular, the multiparticle phase space constraints for a system of 

particles project onto the two-particle space as ll28ll 



Q.{PUP2) =l-Mi- {pi.t-P2.t] -M2 ■ {pi.z-P2.,:} (15) 



where 
Ml 

Mi = 



-M3-{Ei-E2}+M4-{Ei+E2}-^., 



M2 



N{pf) 



N{{E^)-{Ey) 



, M4: 



A^((£2)_(£)2)- 



(16) 



The notation {X} in Equation [Tsl is used to indicate that X 
is the average of a two-particle quantity which depends on pi 
and p2 (or q, etc). In particular, 

T^..^ ^ / d^Pl J d^hP {Pl)P {P2)X5 {q - (pi -P2)) 
^ Jd^piJd^p2P{pi)P{p2)Hq-{pi-p2)) ' 

(17) 

where P represents the single-particle probability first seen in 
Equation[T] 

In practice, this means generating histograms in addition to 
A (q) and B (q) (cf. Equation|2l) as one loops over mixed pairs 
of particles / and j in the data analysis. For example 

{pi,T-P2,T}{q) = ' 

where the sum in the numerator runs over all pairs in all 
events. 

In Equation [15] the four fit parameters M, are directly re- 
lated to five physical quantities, (A^ - the number of particles, 
(Pr)' (P?)' (£■)) through Eq. [16] Assuming that 



{pl) + {p',)+ml 



(19) 
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where is the mass of a typical particle in the system (for 
our pion-dominated system, m^, w mj[), then one may solve for 
the physical parameters. For example, 

N^^l^±^l-^. (20) 

Since we cannot know exactly the values of (E^) etc, that 
characterize the underlying distribution in these collisions, we 
treat the M,- as free parameters in our fits, and then consider 
whether their values are mutually c omp atible and physical. 
For a more complete discussion, see l28Ll47ll . 

In lEsll . the correlations leading to EquationfTslwere called 
"EMCICs" (short for Energy and Momentum Conservation- 
Induced Correlations); we will refer to fits using this function 
with this acronym, in our figures. 

C. Parameter counting 

As mentioned, we will be employing a number of different 
fitting functions, each of which contains several parameters. 
It is appropriate at this point to briefly take stock. 

In essentially all modern HBT analyses, on the order of 
5-6 parameters quantify the femtoscopic correlations. For 
the common Gaussian fit (equation [TTli. one has three "HBT 
radii," the chaoticity parameter, and the normalization 91 . 
Recent "imaging" fits approximate the two-particle emission 
zone as a sum of spline functions, the weights of which are the 
parameters UK]; the number of splines (hence weights) used 
is ^ 5. Other fits (e.g. double Gaussian, exponential-plus- 
Gaussian) ifisi I49I1 contain a similar number of femtoscopic 
parameters. In all cases, a distinct set of parameters is ex- 
tracted for each selection of P (c.f. equation|2]and surrounding 
discussion). 

Accounting for the non-femtoscopic correlations inevitably 
increases the total number of fit parameters. The "^ — p" func- 
tional form (eq.O involves two parameters, the "5 — q" form 
(eq. [TJI i three, and the EMCIC form (eq. [TSl l four However, 
it is important to keep in mind that using the ^ — p (5 — ^) 
form means 2 (3) additional parameters /or each selection of 
P when forming the correlation functions. On the other hand, 
the four EMCIC parameters cannot depend on P. Therefore, 
when fitting (q) for four selections of P, use of the ^ — p, 
5 — ^ and EMCIC forms increases the total number of param- 
eters by 8, 12 and 4, respectively. 

III. ANALYSIS DETAILS 

As mentioned in Section |l] there is significant advantage 
in analyzing p + p collisions in the same way that heavy ion 
collisions are analyzed. Therefore, the results discussed in this 
paper are produced with the same techniques and acceptance 
cuts as have been used for previous pion femtoscopy studies 
by STAR ifsol - lssll . Here we discuss some of the main points; 
full systematic studies of cuts and techniques can be found 
in 1I52I1 . 



The primary sub-detector used in this analysis to recon- 
struct particles is the Time Projection Chamber (TPC) lf54ll . 
Pions could be identified up to a momentum of 800 MeV /c by 
correlating their momentum and specific ionization loss 
(dE/dx) in the TPC gas. A particle was considered to be a 
pion if its dE/dx value for a given momentum was within 
two sigma of the Bichsel expectation ifssll (an improvement 
on the Bethe-Bloch formula [5^ for thin materials) for a pion, 
and more than two sigma from the expectations for electrons, 
kaons and protons. By varying the cuts on energy loss to al- 
low more or less contamination from kaons or electrons, we 
estimate that impurities in the pion sample lead to an uncer- 
tainty in the femtoscopic scale parameters (e.g. HBT radii) 
of only about 1%. Particles were considered for analysis if 
their reconstructed tracks produced hits on at least 10 of the 
45 padrows, and their distance of closest approach (DCA) to 
the primary vertex was less than 3 cm. The lower momentum 
cut of 120 MeV/c is imposed by the TPC acceptance and the 
magnetic field. Only tracks at midrapidity {\y\ < 0.5) were 
included in the femtoscopic analysis. 

Events were recorded based on a coincidence trigger of 
two Beam-Beam Counters (BBCs), annular scintillator detec- 
tors located ±3.5 m from the interaction region and cover- 
ing pseudorapidity range 3.3 < |ri| < 5.0. Events were se- 
lected for analysis if the primary colUsion vertex was within 
30 cm of the center of the TPC. The further requirement that 
events include at least two like-sign pions increases the av- 
erage charged particle multiplicity with |r|| < 0.5 from 3.0 
(without the requirement) to 4.25. Since particle pairs en- 
ter into the correlation function, the effective average mul- 
tiplicity is higher; in particular, the pair-weighted charged- 
particle multiplicity at midrapidity is about 6.0. After event 
cuts, about 5 million minimum bias events from p + p colli- 
sions at y/s ~ 200 GeV were used. 

Two-track effects, such as splitting (one particle recon- 
structed as two tracks) and merging (two particles recon- 
structed as one track) were treated identically as has been done 
in STAR analyses of Au + Au collisions 1531 • Both effects 
can affect the shape of C{q) at very low \q\ < 20 MeV/c, re- 
gardless of the colliding system. However, their effect on the 
extracted sizes in p + p collisions turns out to be smaller than 
statistical errors, due to the fact that small (^ 1 fm) sources 
lead to large (~ 200 MeV/c) femtoscopic structures in the 
correlation function. 

The analysis presented in this paper was done for four bins 
in average transverse momentum kr (= j\ipT,i+PT,2)\)- 
150-250, 250-350, 350-450 and 450-600 MeV/c. The' sys- 
tematic errors on femtoscopic radii due to the fit range, parti- 
cle mis-identification, two-track effects and the Coulomb ra- 
dius (used to calculate K^oui in Eqs.[TOland [TTI) are estimated 
to be about 10%, similar to previous studies li52R . 

IV. RESULTS 

In this section, we present the correlation functions and 
fits to them, using the various functional forms discussed in 
Section The mj- and multiplicity dependence of femto- 
scopic radii from these fits are compared here, and put into 
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the broader context of data from heavy ion and particle colli- 
sions in the next section. 

Figure [T] shows the two-pion correlation function for 
minimum-bias p + p collisions for 0.35 < kj < 0.45 GeV/c. 
The three-dimensional data is represented with the traditional 
one-dimensional Cartesian projections |fTll). For the projec- 
tion on qo, integration in and qi was done over the range 
[0.00.0.12] GeV/c. As discussed in Section HI] and in more 
detail in ll28ll . the full structure of the correlation function is 
best seen in the spherical harmonic decomposition, shown in 
Figs. ED 

In what follows, we discuss systematics of fits to the cor- 
relation function, with particular attention to the femtoscopic 
parameters. It is important to keep in mind that the fits are 
performed on the full three-dimensional correlation function 
C{q). The choice to plot the data and fits as spherical har- 
monic coefficients A/„, or as Cartesian projections along the 
"out," "side" and "long" directions is based on the desire to 
present results in the traditional format (projections) or in a 
representation more sensitive to the three-dimensional struc- 
ture of the data ||28|] . In particular, the data and fits shown in 
Fig lU for fc7-=0. 35-0.45 GeV/c, are the same as those shown 
in Fig. m 

A. Transverse mass dependence of 3D femtoscopic radii 

Femtoscopic scales from three-dimensional correlation 
functions are usually extracted by fitting to the functional form 
given in Equation[TT] In order to make connection to previous 
measurements, we employ the same form and vary the treat- 
ment of non-femtoscopic effects as discussed in Section HTB] 
The fits are shown as curves in Figs. [T]|5j the slightly fluctu- 
ating structure observable in the sensitive spherical harmonic 
representation in Figs. |2]|5]results from finite-binning effects 
in plotting 1,571 . 

Dashed green curves in Figs. [T]|5] represent the "standard 
fit," in which non-femtoscopic correlations are neglected al- 
together (i2 = 1). Black dotted and purple dashed curves, 
respectively, indicate "5 — q" (Equation [T3) and "^ — p" 
(Equation [T4l i forms. Solid red curves represent fits in 
which the non-femtoscopic contributions follow the EMCIC 
(Equation [TSl l form. None of the functional forms perfectly 
fits the experimental correlation function, though the non- 
femtoscopic structure is semi-quantitatively reproduced by the 
ad-hoc 8 — q and ^ — p fits (by construction) and the EMCIC 
fit (non-trivially). Rather than invent yet another ad-hoc func- 
tional form to better fit the data, we will consider the radii 
produced by all of these forms. 

The fit parameters for these four fits, for each of the four 
kj bins, are given in Tables Hi ll VI Considering first the non- 
femtoscopic correlations, we observe that the ad-hoc fit pa- 
rameters do.s.L in TableHJand ^ and p in Table Ulllare different 
for each kj bin. Due to their physical meaning, the EMCIC 
parameters Mi_4 are fixed for all kj values, as indicated in 
Table II VI Setting the characteristic particle mass to that of the 
pion and using Equations [T6l [T9] and l20l the non-femtoscopic 
parameters listed in TablellVlcorrespond to the following val- 




FIG. 1: (Color online) Cartesian projections of the 3D correlation 
function from p + p collisions at ^^=200 GeV for kj = [0.35,0.45] 
GeV/c (blue triangles). Femtoscopic effects are parameterized with 
the form in Eq. [TT] different curves represent various parameteriza- 
tions of non-femtoscopic correlations used in the fit and described in 
detail in Sec HTBl 



ues characteristic of the emitting system: 

A^= 14.3 ±4.7 
ifl) =0.17 ±0.06 (GeV/c)^ 
(pi) =0.32±0.13 (GeV/c)2 
(£2) = 0.51 ±0.11 GeV- 
(E) =0.68 ±0.08 GeV. 

These values are rather reasonable ll47ll . 

HBT radii from the different fits are plotted as a func- 
tion of transverse mass in Fig. |6] The treatment of the non- 
femtoscopic correlations significantly affects the magnitude 
of the femtoscopic length scales extracted from the fit, espe- 
cially in the "out" and "long" directions, for which variations 
up to 50% in magnitude are observed. The dependence of 

the radii on mj = ^kj + rrp- is quite similar in all cases. We 
discuss this dependence further in Section IVl 
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kT [GeV/f] 


Ro [fm] 


R, [fm] 


R, [fm] 


\ 


xVndf 


[0.15,0.25] 
[0.25,0.35] 
[0.35,0.45] 
[0.45,0.60] 


0.84 ±0.02 
0.81 ±0.02 
0.71 ±0.02 
0.68 ±0.02 


0.89±0.01 
0.88±0.01 
0.82±0.02 
0.68±0.01 


1.53 ±0.02 
1.45 ±0.02 
1.31±0.02 
1.05 ±0.02 


0.422 ±0.004 
0.422 ±0.005 
0.433 ±0.007 
0.515±0.009 


2012/85 
1852/85 
941/85 
278 / 85 



TABLE I: Fit results from a fit to data from p-\-p collisions at y^= 200 GeV using Eq. [TT]to parameterize the femtoscopic correlations 
("standard fit"). 



kT [GeV/c] 


Ro [fm] 


R, [fm] 


R, [fm] 


\ 


5„ 


8., 


s, 


xVndf 


[0.15,0.25] 
[0.25,0.35] 
[0.35,0.45] 
[0.45,0.60] 


1.30±0.03 
1.21 ±0.03 
1.10±0.03 
0.93 ±0.03 


1.05 ±0.03 
1.05 ±0.03 
0.94 ±0.03 
0.82 ±0.03 


1.92 ±0.05 
1.67 ±0.05 
1.37 ±0.05 
1.17 ±0.05 


0.295 ±0.004 
0.381 ±0.005 
0.433 ±0.007 
0.480 ±0.009 


0.0027 ±0.0026 
0.0201 ±0.0054 
0.0457 ±0.0059 
0.0404 ±0.0085 


-0.1673 ±0.0052 
-0.1422 ±0.0051 
-0.0902 ±0.0053 
-0.0476 ± 0.0093 


-0.2327 ±0.0078 
-0.2949 ±0.0081 
-0.2273 ±0.0090 
-0.1469 ±0.0104 


471/82 
261/82 
251/82 
189/82 



TABLE IL Fit results from a fit to data from p-\- p collisions at 200 GeV using Eq.[TT]to parameterize the femtoscopic correlations and 
Eq.[T3]for non-femtoscopic ones ("8-5 fit"). 
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FIG. 2: (Color online) The first three non-vanishing moments of the 
spherical harmonic decomposition of the correlation function from 
p + p collisions at VJ=200 GeV, for kj = [0.15,0.25] GeV/c. Fem- 
toscopic effects are parameterized with the form in Eq. 1111 differ- 
ent curves represent various parameterizations of non-femtoscopic 
correlations used in the fit and described in detail in Sec. IIIBI The 
superscript "RE" in Af^^^ stands for the real part of A; ,„ . 



B. Transverse mass and multiplicity dependence of ID 
femtoscopic radii 



Since three-dimensional correlation functions encode more 
information about the homogeneity region than do one- 
dimensional correlation functions, they are also more statis- 
tics hungry. Therefore, most previous particle physics experi- 
ments have constructed and analyzed the latter For the sake of 
making the connection between our results and existing world 
systematics, we perform similar analyses as those found in the 
literature. 

The first important connection to make is for the nij- 




A CF 

standard fit 

h-q fit 
C-Pfit 
EMCIC fit 



0.2 ^ 0.4 
\q\ [GeV/c] 



FIG. 3: (Color online) As for Fig.lS but for kj = [0.25, 0.35] GeV/c 



dependence of HBT radii from minimum-bias /? + p colli- 
sions. We extract the one-dimensional HBT radius R^^ asso- 
ciated with the femtoscopic form in Equation [TOl using three 
forms for the non-femtoscopic terms. For four selections in 
kj. Table W\ lists the fit parameters for the "standard" fit that 
neglects non-femtoscopic correlations altogether (il = 1). Ta- 
bles|VI|and |VII| list results when using the 1 -dimensional 8 — ^ 
form (Equation [T2b and the EMCIC form (Equation [TSll, re- 
spectively. In performing the EMCIC fit, the non-femtoscopic 
parameters Mi_4 were kept fixed at the values listed in Ta- 
ble |IVl 

The one-dimensional radii from the three different treat- 
ments of non-femtoscopic effects are plotted as a function of 
niT in Fig. [T] The magnitude of the radius using the ad-hoc 
8 — ^ fit is ^25% larger than that from either the standard or 
EMCIC fit, but again all show similar dependence on mj-. 

In order to compare with the multiplicity dependence of 
^7- -integrated HBT radii reported in high energy particle col- 
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kj [GeV/c] 


Ro [fm] 


R, [fm] 


R, [fm] 


X 




P 


xVndf 


[0.15,0.25] 
[0.25,0.35] 
[0.35,0.45] 
[0.45,0.60] 


1.24 ±0.04 
1.14 ±0.05 
1.02 ±0.04 
0.89 ±0.04 


0.92±0.03 
0.89 ±0.04 
0.81±0.05 
0.71 ±0.05 


1.71±0.04 
1.37 ±0.08 
1.20 ±0.07 
1.09 ±0.06 


0.392 ±0.008 
0.378 ±0.006 
0.434 ±0.008 
0.492 ±0.009 


0.0169±0.0021 
0.0193 ±0.0034 
0.0178 ±0.0029 
0.01 14± 0.0023 


-0.0113±0.0019 
-0.0284 ±0.0031 
-0.0289 ±0.0032 
-0.0301 ±0.0041 


1720 / 83 
823 / 83 
313/83 
190 / 83 



TABLE III: Fit results from a fit to data from p + p collisions at y(s= 200 GeV using Eq.[TT]to parameterize the femtoscopic correlations and 
Eq.[T4]for non-femtoscopic ones ('% — [3 fit"). 



kj [GeV/c] 


Ro [fm] 


Rs [fm] 


R, [fm] 


X 


Ml (GeV/c)-2 


Mo (GeV/c)"2 


M-i GeV-2 


M4 GeV-' 


X-/ndf 


[0.15,0.25] 
[0.25,0.35] 
[0.35,0.45] 
[0.45,0.60] 


1.06±0.03 
0.96 ±0.02 
0.89 ±0.02 
0.78 ±0.04 


1.00 ±0.04 
0.95 ±0.03 
0.88 ±0.02 
0.79 ±0.02 


1.38 ±0.05 
1.21 ±0.03 
1.08 ±0.04 
0.94 ±0.03 


0.665 ±0.005 
0.588 ±0.006 
0.579 ±0.009 
0.671 ±0.028 


0.43 ±0.07 


0.22 ±0.06 


1.51±0.12 


1.02 ±0.09 


2218/336 



TABLE IV: Fit results from a fit to data from p-\- p collisions at v^= 200 GeV using Eq.[TT]to pai^ameterize the femtoscopic correlations and 
Eq.[T5]for non-femtoscopic ones ("EMCIC fit"). 




kj [GeV/c] 


R\m [fm] 


X 


xVndf 


[0.15,0.25] 
[0.25,0.35] 
[0.35,0.45] 
[0.45,0.60] 


1.32±0.02 
1.26 ±0.02 
1.18±0.02 
1.05 ±0.03 


0.345 ±0.005 
0.357 ±0.007 
0.348 ±0.008 
0.413±0.012 


265 / 27 
203 / 27 
243 / 27 
222 / 27 



TABLE V: Fit results from a fit to ID correlation function from 
p-\- p collisions at -^J = 200 GeV using Eq. |6] to parameterize the 
femtoscopic correlations ("standard fit"). 



kj [GeV/c] 


^inv [fm] 


X 


S 


xVndf 


[0.15,0.25] 
[0.25,0.35] 
[0.35,0.45] 
[0.45,0.60] 


1.72 ±0.04 
1.65 ±0.04 
1.49 ±0.05 
1.41 ±0.06 


0.285 ±0.007 
0.339 ±0.009 
0.308±0.011 
0.338±0.016 


0.237 ±0.007 
0.163 ±0.008 
0.180±0.015 
0.228 ±0.017 


86/26 
80/26 
71/26 
78/26 



TABLE VI: Fit results from a fit to ID correlation function from 
p + p collisions at = 200 GeV using Eq. [6] to parameterize the 
femtoscopic correlations and Eq.[T2]for non-femtoscopic ones ("S — 
q fit"). 



lisions, we combine fcj- bins and separately analyze low 
(dNch/d\\ < 6) and high {dNch/d\\ > 7) multiplicity events. 
The choice of the cut was dictated by the requirement of suf- 
ficient pair statistics in the two event classes. Fit parame- 
ters for common fitting functions are given in Table IVIIII for 
minimum-bias and multiplicity-selected collisions. 

Figure[8]shows the multiplicity dependence of the common 



one-dimensional HBT radius R[^^, extracted by parameteriz- 
ing the femtoscopic coiTelations according to Equation [TO] 
Non-femtoscopic effects were either ignored ("standard fit" 
n = 1) or parameterized with the "5 — q" (Eq. [T2b or EM- 
CIC (Eq. [TSl l functional form. In order to keep the parameter 
count down, the EMCIC kinematic parameters {{pj), {p^). 
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FIG. 6: (Color online) The raj-dependence of the 3D femtoscopic 
radii in p + p collisions at ^/s = 200 GeV for different parameteriza- 
tions of the non-femtoscopic correlations. See text for more details. 
Data have been shifted slightly in the abscissa, for clarity. 
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FIG. 7: (Color online) The otj- -dependence of ^inv from p + p col- 
lisions at = 200 GeV for different parameterizations of the non- 
femtoscopic correlations used in the fit procedure. 



kj [GeV/c] 


^inv [fm] 


X 


X-/ndf 


[0.15,0.25] 
[0.25,0.35] 
[0.35,0.45] 
[0.45,0.60] 


1.38±0.03 
1.32±0.03 
1.23 ±0.04 
1.14±0.05 


0.347 ±0.005 
0.354± 0.006 
0.349 ±0.009 
0.411±0.013 


99/27 
97/27 
86/27 
80/27 



TABLE VII: Fit results from a fit to ID correlation function from 
p + p collisions at ^/s = 200 GeV using Eq. [6] to parameterize the 
femtoscopic correlations and Eq.[T5]for non-femtoscopic ones ("EM- 
CIC fit"). The non-femtoscopic parameters Mi_4 were not varied, 
but kept fixed to the values in Table II V I 



scales, both with kj and multiplicity, are robust. 

Table llXl lists fit parameters to two-dimensional correlation 
functions in qj and qo, using Equations [7] and [8] The radius 
from the former fit is approximately twice that of the latter, as 
expected (c.f. Sec. Ill Al ). These values will be compared with 
previously measured data in the next section. 



I 1.5 

> 
c 

^~ 1 
0.5 



















• 

















• 




★ 








it 


■ 






* 


□ 








■ 
















min-bias 








■ Standard fit 


□ 








• b-q fit 


o 








★ EMCIC fit 





10 



FIG. 8: (Color online) The multiplicity dependence of Ri„y from 
p + p collisions at = 200 GeV for different parameterizations of 
the non-femtoscopic correlations. Pions within the range of kj = 
[0.15,0.60] GeV/c were used in the analysis. 



V. COMPARISON WITH WORLD SYSTEMATICS 



(E^), {£)) were kept fixed to the values obtained from the 3- 
dimensional fit, and only N was allowed to vary. In all cases, 
Rinv is observed to increase with multiplicity. Parameterizing 
non-femtoscopic effects according to the EMCIC form gives 
similar results as a "standard" fit ignoring them, whereas the 
"6 — q" form generates an offset of approximately 0.3 fm, sim- 
ilar to all three- and one-dimensional fits discussed above. 
That different numerical values are obtained for somewhat 
different fitting functions, is not surprising. The point we fo- 
cus on is that the systematic dependences of the femtoscopic 



In this section, we make the connection between femto- 
scopic measurements in heavy ion collisions and those in par- 
ticle physics, by placing our results in the context of world 
systematics from each. 

A. Results in ttie Context of Heavy Ion Systematics 

The present measurements represent the first opportunity to 
study femtoscopic correlations from hadronic collisions and 
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method 


fit parameter 


{dN,iJdr\} 


4.25 (min-bias) 


3.47 


8.75 




^inv [fm] 


1.21±0.01 


1.09 ±0.02 


1.34±0.02 


standard fit 


X 


0.353 ±0.003 


0.347 ±0.04 


0.356±0.03 




xVndf 


202 / 27 


100 / 27 


92/27 




^inv [fm] 


1.61±0.01 


1.50±0.03 


1.76±0.03 


6-^ fit 


X 


0.312±0.003 


0.275 ±0.005 


0.322 ±0.007 




§Ginv [c/GeV] 


-0.191 ±0.003 


-0.242 ±0.005 


-0.194±0.006 




xVndf 


159/26 


83/26 


73/26 




^inv [fm] 


1.32±0.02 


1.22±0.03 


1.46 ±0.02 


EMCIC fit 




0.481 ±0.003 


0.485 ±0.003 


0.504 ±0.004 




N 


14.3±4.7 


11.8±7.1 


26.3 ±8.4 




xVndf 


161/26 


80/26 


75/26 



TABLE VIII: Multiplicity dependence of fit results to ID correlation function from p + p collisions at ^/s — 200 GeV for different fit parame- 
terizations. 



method 


fit parameter 


{dN,iJdr]) 


4.25 (min-bias) 


3.47 


8.75 




Rb [fm] 


1.79±0.01 


1.61±0.02 


1.92 ±0.02 


Eq.m 


T [fm/c] 


1.03 ±0.02 


0.98 ±0.02 


1.24 ±0.03 




I 


0.353 ±0.003 


0.354± 0.003 


0.334± 0.004 




xVndf 


5308 / 896 


2852 / 896 


1890/896 




Rg [fm] 


1.0I±0.01 


0.89±0.0I 


1.07±0.01 


Eq.l 


T [fm/c] 


0.76±0.01 


0.73 ±0.02 


0.91 ±0.02 


X 


0.353 ±0.003 


0.352 ±0.003 


0.332 ±0.004 




xVndf 


5749 / 896 


3040 / 896 


2476 / 896 



TABLE IX: Multiplicity dependence of fit parameters to two-dimensional correlation functions from p + p collisions at ^/s = 200 GeV using 
Equations|7]and[8l To consistently compare to previous measurements, £1 was set to unity (c.f. Equation[5]l. 



heavy ion collisions, using the same detector, reconstruction, 
analysis and fitting techniques. The comparison should be di- 
rect, and differences in the extracted HBT radii should arise 
from differences in the source geometry itself. In fact, espe- 
cially in recent years, the heavy ion community has generally 
arrived at a consensus among the different experiments, as far 
as analysis techniques, fitting functions and reference frames 
to use. This, together with good documentation of event se- 
lection and acceptance cuts, has led to a quantitatively consis- 
tent world systematics of femtoscopic measurements in heavy 
ion collisions over two orders of magnitude in collision en- 
ergy iflltl ; indeed, at RHIC, the agreement in HBT radii from 
the different experiments is remarkably good. Thus, inas- 
much as star's measurement of HBT radii from p + p colli- 
sions may be directly compared with STAR'S HBT radii from 
Au + Au collisions, they may be equally well compared to the 
world's systematics of all heavy ion collisions. 

As with most heavy ion observables at low transverse mo- 
mentum ifssll . the HBT radii Rs and Ri scale primarily with 
event multiplicity ifTTll (or, at lower energies, with the number 
of particles of different species lf59ll60ll ) rather than with en- 
ergy or impact parameter. The radius which nontrivially 
combines space and time, shows a less clear scaling ifTHl . re- 



taining some energy dependence. As seen in Fig. |9] the radii 
from p + p collisions at y/s ^ 200 GeV fall naturally in line 
with this multiplicity scaling. On the scale relevant for this 
comparison, the specific treatment of non-femtoscopic coiTe- 
lations is unimportant. 

One of the most important systematics in heavy ion fem- 
toscopy is the mj-dependence of HBT radii, which di- 
rectly measures space-momentum correlations in the emit- 
ting source at freeze-out; in these large systems, the mj- 
dependence is often attributed to collective flow f^j. As 
we saw in Fig. |6] a significant dependence is seen also for 
p + p collisions. Several authors [e.g. [H [M El [H, [Ml 
have remarked on the qualitative "similarity" of the mr- 
dependence of HBT radii measured in high energy particle 
collisions, but the first direct comparison is shown in Fig. [TO] 
There, the ratios of the three dimensional radii in Au + Au col- 
lisions to p + p radii obtained with different treatments of the 
non-femtoscopic correlations, are plotted versus mj-. Well be- 
yond qualitative similarity, the ratios are remarkably flat- i.e. 
the nij- -dependence mp + p collisions is quantitatively almost 
identical to that in Au + Au collisions at RHIC. We speculate 
on the possible meaning of this in Section lVB] 
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B. Results in the context of high-energy particle measurements 
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FIG. 9: (Color online) The multiplicity dependence of the HBT radii 
from p + p , Cu + Cu fssll and Au + Au (s^ Issll collisions from 
STAR compared with results from other experiments ITlll . Left and 
right panels show radii measured with (kj) ~ 0.2 and 0.39 GeV/c, 
respectively. Radii from p + p collisions are shown by blue ("stan- 
dard fit") and red ("EMCIC fit") stars. 



• Au+Au {0-5%)/ p+p [St] 
T Au+Au (0-5%) /p+p [6-?] 
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A Au+Au (50-80%) / p+p l&-q] 

□ Au+Au (50-80%) / p+p [C-p] 
Au+Au (50-80%) /p+p [EMCIC] 
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FIG. 10: (Color online) The ratio of the HBT radii from Au + Au col- 
lisions |52|| to results from p + p collisions plotted versus the trans- 



verse mass. 



Recently, a review of the femtoscopic results 112011 from par- 
ticle collisions like p + p , p+ p and e+ + studied at differ- 
ent energies has been published. Here, we compare STAR re- 
sults from p + p collisions at ^/s — 200 GeV with world sys- 
tematics. 



STARp+p @ 200 GeV 

■ standard fit 

★ EMCIC fit 

■ E735P+P @1.8TeV 

ABCDHW 

• p+p @ 44 GeV 

O p+p @ 62 GeV 




5 10 15 20 



FIG. II: (Color online) The multiplicity dependence of the ID fem- 
toscopic radius i^inv from hadronic collisions measured by STAR, 
E735 m, and ABCDHW lH collaborations. 

The multiplicity dependence of femtoscopic parameters 
from one- and two-dimensional correlation functions are 
shown in Figs. [TTI and [T2l For any given experiment, the ra- 
dius parameter increases with event multiplicity. However, 
in contrast to the nearly "universal" multiplicity dependence 
seen in heavy ion collisions (c.f. Fig. |9]l, only a qualitative 
trend is observed, when the different measurements are com- 
pared. 

There are several possible reasons for this lack of "univer- 
sality" lEoll . Clearly one possibility is that there is no universal 
multiplicity dependence of the femtoscopic scales; the under- 
lying physics driving the space-time freezeout geometry may 
be quite different, considering y/s varies from 44 to 1800 GeV 
in the plot. However, even if there were an underlying univer- 
sality between these systems, it is not at all clear that it would 
appear in this figure, due to various difficulties in tabulating 
historical data 112011 . Firstly, as discussed in Section HIl the ex- 
periments used different fitting functions to extract the HBT 
radii, making direct comparison between them difficult. Sec- 
ondly, as we have shown, the radii depend on both multiplic- 
ity and kr- Since, for statistical reasons, the results in Fig.|9] 
are integrated over the acceptance of each experiment, and 
these acceptances differ strongly, any universal scaling would 
be obscured. For example, since the acceptance of Tevatron 
experiment E735 jSfil] is weighted towards higher kj than the 
other measurements, one expects a systematically lower HBT 
radius, at a given multiplicity. Indeed, even the "universal" 
multiplicity scaling in heavy ion collisions is only universal 
for a fixed selection in kj- Thirdly, the measure used to quan- 
tify the event multiplicity varies significantly in the historical 
literature; thus the determination of {dNch/dr[) for any given 
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• AFSp+p@53GeV(ii*) 

• AFSp+p (3) 53-63 GeV (ti*) 

• AFS a+a @ 1 26 GeV (n*) 
^STARp+p @ 200 GeV 

□ NA5 p+p- @ 200 GeV 

■ NA5 p+p <a> 200 GeV 
A UA1p+p@630GeV 

«G 

if STAR p+p @ 200 GeV 

• UAIp+p @200GeV 
A UA1p+^@630GeV 
T UA1p+^@900GeV 

■ E735P+P @1.8TeV 



FIG. 12: (Color online) The multiplicity dependence of radius and 
timescale parameters to 2-dimensional correlation functions mea- 
sured by STAR, E735 (H, UAl (H, AFS (H and NA5 jH. The 
legend on the right indicates that the first 7 sets of datapoints come 
from fits to Eq.|7] in which case the parameter Rb/2 is plotted in the 
upper panel; the last 5 sets of datapoints come from fits to Eq.[8] for 
which Re is plotted. As discussed in Section Hi Al and confirmed by 
STAR and UAI, Rc ^ Rb/2. The UAl Collaboration set x = in 
their fits. 



experiment shown in Fig.|9]is only approximate. 

From the discussion above, we cannot conclude definitively 
that there is- or is not- a universal multiplicity scaling of fem- 
toscopic radii in high energy hadron-hadron collisions. We 
conclude only that an increase of these radii with multiplicity 
is observed in all measurements for which y/s > 40 GeV and 
that the present analysis of p + p collisions is consistent with 
world systematics. 

In Section IIVI we discussed the pr-dependence of HBT 
radii observed in our analysis. Previous experiments on 
high-energy collisions between hadrons- and even leptons- 
have reported similar trends. As discussed above, direct 
comparisons with historical high-energy measurements are 
problematic. Comparisons between fit parameters to 1- and 
2-dimensional correlation functions are shown in Figs. [13] 
and [141 All experiments observe a decrease in femtoscopic 
parameters with increasing transverse momentum. Our radii 
at y/s=2QQ GeV fall off similarly or somewhat more than 
those measured at an order of magnitude lower energy at the 
SPS ll30[|43ll . and less than those measured at an order of mag- 
nitude higher energy at the Tevatron ll36ll . It is tempting to 
infer that this compilation indicates an energy evolution of the 
/?r -dependence of femtoscopic radii. However, given our pre- 
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FIG. 13: (Color online) One-dimensional femtoscopic radii from 
p + p collisions at RHIC and p + p collisions at the Tevatron (3^ 
are plotted versus the transverse momentum Pj = [pij + pij) /2 
(c.f. Eq.H. 
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FIG. 14: (Color online) The transverse momentum dependence of fit 
parameters to two-dimensional correlation functions. STAR results 



from fit to Equation [8 
NA27 (H and NA22 



compared to measurements by E735 ya], 
I el). The SPS experiments NA22 and NA27 
set T = in their fits. STAR and E735 data plotted versus Pj = 
{pl.T + P2.t) /2 (c.f. Eq.O. NA27 reported results in terms of \P\ 
and NA22 in terms of 2\P\. For purposes of plotting here, Pj = 
•\/(2/3)|F| was assumed. 



vious discussion, we conclude only that there is qualitative 
agreement between experiments at vastly different collision 
energies, and all show similar pr dependence. 

Systematics in 3-dimensional HBT radii from hadron col- 
lisions are less clear and less abundant, though our measure- 
ments are again qualitatively similar to those reported at the 
SPS, as shown in Fig. [15] There, we also plot recent results 
from e+ + collisions at LEP; in those 3-dimensional analy- 
ses, the "longitudinal" direction is the thrust axis, whereas the 
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beam axis is used in hadron-hadron collisions, as in heavy ion 
collisions. 



STARp+p @ 200 GeV (stand.) 
STARp+p @ 200 GeV (EMCICs) 
NA22 ;i+p@21.7GeV 
NA49* p+p @ 17.3 GeV (prelim) 
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L3* e*+c- @ 91 GeV 
DELPHI e*+e- @ 91 GeV 
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FIG. 15: (Color online) The transverse mass dependence of 3D fem- 
toscopic radii from particle collisions. Data from NA22 (soll . NA49 
preliminary lH, OPAL (H, L3 (H, DELPHI (H. 



VI. DISCUSSION 

We have seen that HBT radii from p + p collisions at RHIC 
are qualitatively consistent with the trends observed in parti- 
cle collisions over a variety of collision energies. Further, they 
fall quantitatively into the much better-defined world system- 
atics for heavy ion collisions at RHIC and similar energies. 
Particularly intriguing is the nearly identical dependence on 
niT of the HBT radii in p + p and heavy ion collisions, as this 
dependence is supposed ll24ll69ll to reflect the underlying dy- 
namics of the latter Several possible sources of an nij- depen- 
dence of HBT radii in small systems have been put forward to 
explain previous measurements. 

1. Alexander ef al. iItoIItiIi have suggested that the Heisen- 
berg uncertainty principle can produce the transverse momen- 
tum dependence of femtoscopic radii in e+ + e^ collisions. 
However, as discussed in 1I20I1 . a more detailed study of the 
results from e+ + collisions complicates the quantitative 
comparisons of the data from various experiments and thus 
the interpretation. Additionally, the arguments from iFtoI FtHi 
apply only to the longitudinal direction (/?/), so could not ex- 
plain the dependence of all three radii. 

2. In principle, string fragmentation should also gener- 
ate space-momentum correlations in small systems, hence an 
nij dependence of the HBT radii. However, there are almost 
no quantitative predictions that can be compared with data. 
The numerical implementation PYTHIA, which incorporates 



the Lund string model into the soft sector dynamics, imple- 
ments Bose-Einstein enhancement only as a crude parameter- 
ization designed to mock up the effect [c.f. Section 12.4.3 of 
P72I1 for the purpose of estimating distortions to W-boson in- 
variant mass spectrum. Any Bose-Einstein correlation func- 
tion may be dialed into the model, with 13 parameters to 
set the HBT radius, lambda parameter, and correlation shape; 
there is no first-principles predictive power On more general 
grounds, the mass dependence of the femtoscopic radii cannot 
be explained within a Lund string model ll73l - l75ll . 

3. Long-lived resonances may also generate the space- 
momentum dependence of femtoscopic radii ll76ll . How- 
ever, as discussed in | |20|) . the resonances would affect the 
HBT radii from p + p collisions differently than those from 
Au + Au collisions, since the scale of the resonance "halo" is 
fixed by resonance lifetimes while the scale of the "core" is 
different for the two cases. Thus it would have to be a co- 
incidence that the same mr dependence is observed in both 
systems. Nevertheless, this avenue should be explored further 

4. Bialas et al. have introduced a model iItsIi based on a di- 
rect proportionality between the four-momentum and space- 
time freeze-out position; this model successfully described 
data from e+ + collisions. The physical scenario is based 
on freezeout of particles emitted from a common tube, af- 
ter a fixed time of 1.5 fm/c. With a very similar model, 
Humanic iFttIi was able to reproduce femtoscopic radii mea- 
sured at the Tevatron lf36ll only with strong additional hadronic 
rescattering effects. With rescattering in the final state, both 
the multiplicity- and the mj- -dependence of the radii were re- 
produced IttIi . 

5. It has been suggested H [M [H [M III that the pr- 
dependence of HBT radii in very small systems might reflect 
bulk collective flow, as it is believed to do in heavy ion colli- 
sions. This is the only explanation that would automatically 
account for the nearly identical /jj- -scaling discussed in Sec- 
tion [VA] However, it is widely believed that the system cre- 
ated in p + p collisions is too small to generate bulk flow. 

The remarkable similarity between the femtoscopic system- 
atics in heavy ion and hadron collisions may well be coinci- 
dental. Given the importance of the mj-dependence of HBT 
radii in heavy ion collisions, and the unclear origin of this 
dependence in hadron collisions, further theoretical investiga- 
tion is clearly called for. Additional comparative studies of 
other soft-sector observables (e.g. spectra) may shed further 
light onto this coincidence. 



VII. SUMMARY 

We have presented a systematic femtoscopic analysis of 
two-pion correlation functions from p + p collisions at RHIC. 
In addition to femtoscopic effects, the data show correlations 
due to energy and momentum conservation. Such effects have 
been observed previously in low-multiplicity measurements at 
Tevatron, SPS, and elsewhere. In order to compare to histor- 
ical data and to identify systematic effects on the HBT radii, 
we have treated these effects with a variety of empirical and 
physically-motivated formulations. While the overall magni- 
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tude of the geometric scales vary with the method, the impor- 
tant systematics do not. 

In particular, we observe a significant positive correlation 
between the one- and three-dimensional radii and the multi- 
plicity of the collision, while the radii decrease with increas- 
ing transverse momentum. Qualitatively, similar multiplicity 
and momentum systematics have been observed previously in 
measurements of hadron and electron collisions at the SppS, 
Tevatron, ISR and LER However, the results from these ex- 
periments could not be directly compared to those from heavy 
ion collisions, due to differences in techniques, fitting meth- 
ods, and acceptance. 

Thus, the results presented here provide a unique possibility 
for a direct comparison of femtoscopy in p + p and A + A col- 
lisions. We have seen very similar pr and multiplicity scaling 
of the femtoscopic scales in /9 + /5 as in A + A collisions, inde- 
pendent of the fitting method employed. Given the importance 
of femtoscopic systematics in understanding the bulk sector 
in Au + Au collisions, further exploration of the physics be- 
hind the same scalings in p + p collisions is clearly important, 
to understand our "reference" system. The similarities ob- 
served could indicate a deep connection between the underly- 
ing physics of systems with size on order of the confinement 
scale, and of systems much larger Similar comparisons will 



be possible at the Large Hadron Collider, where the higher 
collision energies will render conservation laws less impor- 
tant, especially for selections on the very highest-multiplicity 
collisions. 
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